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Abstract 

We consider the combined Yang Mills-Dilaton-Gravity system in the presence of 
a Gauss-Bonnet term as it appears in the 4D Effective Superstring Action. We 
give analytical arguments and demonstrate numerically the existence of black hole 
solutions with non-trivial dilaton and Yang Mills hair for the particular case of SU(2) 
gauge fields. The thermodynamical properties of the solutions are also discussed. 
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1. Dilaton gravity with a Gauss-Bonnet term and an SU(2) Yang Mills 
potential. 



The effective theory of gravity resulting from String Theory [Q] at low energies 
includes important modifications of Einstein's theory due to the presence of the 
extra degrees of freedom such as dilatons, axions and Yang Mills fields. The loop- 
corrected Superstring Effective Action through the contribution of the Gauss-Bonnet 
term leads to the existence of singularity-free cosmological solutions [Q || as well as 
to the existence of new dilatonic black holes [J5J. The black hole solutions found 
possess hair outside the horizon as anticipated from general and perturbative in 
a' considerations M IBfl- The existence of black holes with non-trivial hair has been 



demonstrated in the Einstein-Maxwell-Dilaton system || ||T0[ and the Einstein- Yang 
Mills(-Dilaton-Axion) system [pT] . In the present article we shall extend the analysis 
of reference Q by including an SU(2) Yang Mills field. We find new black hole 
solutions of the combined gravity-dilaton-SU(2) Yang Mills system in the (crucial) 
presence of the Gauss-Bonnet term. 



Let us consider the heterotic string effective action, ignoring for simplicity moduli 
and axion fields. Following the notation of reference [|J, we have 




where a' is the Regge slope, g 2 is some gauge coupling constant and F^ u is the Yang 
Mills field strength. The Gauss- Bonnet term is defined as 

Rgb = R^v P aR^ vpa — 4R IMU R' 1U + R 2 (2) 

We shall consider only the contribution of an SU(2) Yang Mills field, assuming trivial 
values for all other Yang Mills as well as "matter" fields. 



At this point we shall make a spherically symmetric ansatz for the space-time 
metric 

ds 2 = -e r dt 2 + e A dr 2 + r 2 (d6 2 + sin 2 6dip 2 ) (3) 
The functions T and A depend solely on the radius r. For the SU(2) Yang Mills 



potential we shall adopt the magnetic ansatz [JTT 

A = 1 + W [—f v dd + Tgsinddip] (4) 

where e is an SU(2) gauge coupling constant, w = w{r) is a radial function, and 
(f r , fg, fip) is the ant i- Hermit ian SU(2) basis expressed in polar coordinates ; e.g. 
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f r — r • r, and [r a ,iV] = e a ft c f c with the indices ranging over (r,8,(p). The field 
strength corresponding to is 

w' (1 — w 2 ) 
F = — [—ftp dr Add + fg sm0 dr A dip] - f r sinO dd A dip (5) 

2. Some analytical considerations. 

Using the static, spherically symmetric ansatz (^|) and (|5]), the dilaton and Yang- 
Mills equation as well as the (tt), (rr) and (88) component of the Einstein's equations 
resulting from ([!]) take the form 



T' - A' 2 n'e^ ( V 

*" + *•(— + = ^(ri' e -' + (i- t -")[r +7 (r-A')i 



^ +w( — 5 — + 0) + - A -^ J -e A = (7) 



o/e 



•V ( 1 + ^(1 - 3e" A ) I = ^ + ^ + + 



q ; e^ /V (1-w 2 ) 2 
2g 2 r V e 2 2e 2 r 2 



L(r'-A') + - - = -^ + — (■•"•" ■ •■" •■ /2 >'-' 

2 r 2 g 2 r 

+ ^(r'-3A')U™(l-- 2 ) 2 e A (10) 



2g 2 e 2 r 4 



The corresponding components of T„ v are 
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Due to the presence of the higher curvature terms, the assumption of positive def- 
initeness of the time- component of the "energy-momentum" tensor breaks down. 
This point is crucial for the evasion of the no-hair conjecture || and the existence 
of new black hole solutions. 



Let us study now our system at infinity. Demanding asymptotically flat solutions, 
we expand the functions e A( * r \ e r (' r \ <p(r) and w(r) in a power series in 1/r and 
substitute them in the equations (^) - fllPl). From eq.(^) in the first order we take 
the constraint: (1 — w^) =0. As a result, the only two options allowed are 
either = or = ±1. The first leads to W\ = in the second order, to W2 = 
in the third order and eventually to the vanishing of all u> n 's. This case corresponds 
to an Abelian potential 

A^^B r = F 9v ^ 1 J± (12) 

with magnetic charge 1/e. The other option, namely = ±1, leads to non- 
vanishing WnS and corresponds to the non- Abelian sector of SU(2). In this case 
the radial magnetic field is B r ~ 1/r 3 implying that the Yang Mills charge Qym 
vanishes at infinity. We shall study the non- Abelian case in this article and, thus, 
make the choice Woo = ±1 which leads us to the following asymptotic form of the 
solution near infinity 

eAW = 1 + 2M + 16M ! ^ + 0( l ) 



e r « = 1- — + 0(1) (14) 

D DM n , 1 . 
4>{r) = 0OO + - + — + 0{-) 15 

w(r) = ±(l + ^i + 0(i)) (16) 
M is the ADM mass and D is the dilaton charge defined over a two-sphere at infinity 

D = -^J d 2 ^V^ (17) 



as 



Let us now direct our attention to the other end of the allowed spatial range, i.e. 
near the event horizon. We note that equation (|9]) can be solved for e A as 

e = Ya (18) 
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where 5 = ±1 and 



a'e* (1-w 2 ) 2 3a'e 4 



ol = 1 — , *y = r 

Ag 2 r 2 e 2 2g 2 

d' 2 r 2 a'e^d' a'e 4 * w' 2 

3 = — — - i -r'(r + ——) + —— — (19) 

P 4 ^ r + 2 g 2 } + 2g 2 e 2 1 j 

As a result e A , as well as A', can be eliminated from equations (^j), (|7|), (|S|) and (|TT|). 
Choosing any three of them, since only three of them are independent, we obtain 
the system 

<t>" =~ (20) 

r" = ~ (21) 

a 

w" (22) 



where d, di, c^, d 3 are complicated functions of 4>', T', w', 0, w, and r. Demanding 
that 4>'l and w'^ are finite, we obtain from ( p0|) and ( p2|) the constraints 

<^ 2 [ 2e 2 ^C h r h (e 2 ^ + e^r 2 + r£) - 8e^r 7 h } + 

<j>' h [e^C 2 + 4e^C h r 2 (e 2 ^ + 2eV 2 + r A h ) - 16rg ] - (23) 

e^ h Clr h + 24e 2 ^C /l r}* + 8e* h C fc r£ - 48e* fc rj| = 

Wfe (l- W 2 )(l + eN4/2r fe ) 

r* (1 - e^C,/4r 2 ) [M) 

where Ch — (1 — w 2 ) 2 /e 2 . Since a'/g 2 always multiplies e^, we may eliminate it 
from our calculations and restore it in the end. Substituting these two into eq.(pl|) 
gives 

r" = -r' 2 + 0m r' = — *— + o(i) (25) 

r - r h 

Expanding equation (|18|) near the horizon in powers of T', for 5 = +1 we obtain 



and 



e 



A 



2r 2 (2r + efy') 
4r 2 _ 



r' + 0(l) (26) 



while the choice 5 = — 1 leads to e A = 0(1) which is not a black hole solution. 
Taking into account all the above, we may conclude that the unique black hole 
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solution with <ft', <ft", w', w" finite and r' — ► oo near the horizon has the expansion 
e r(r) = li(r-r h ) + 0(r-r h ) 2 (27) 
e" A M = \^ r - rh) + 0{r-r h ) 2 (28) 

0(r) = 4> h + (P' h (r - r h ) + ^ (r - r fe ) 2 + C(r - r h f (29) 
w(r) = w h + < (r - r h ) + -y (r - r h f + £>(r - r h f (30) 



where r^, 0^, Wh and 71 are free parameters while from eq. (|26|) we get 



4r 2 _ e 4> h Ch 

Al = 2r 2 {2r h + e*»V h ) (31) 



Note that the constraint fl2"3"|) is an algebraic second order equation for (j>' h which 
has two real solutions, 0' ± , only if its discriminant is positive. This ultimately gives 

(C h -z+)(C h -z-) >0 (32) 

where 



z± = — [-(2 + 3x + 3x 2 ) ± J (2 + 3x + 3x 2 ) 2 + Qx 2 - 1 ] (33) 
or v 

having set x = e^ h j g 2 r\. When x < l/Vo, or 



< ^ (34) 



the discriminant is always positive. This has the obvious interpretation that if the 
Gauss-Bonnet effective coupling a'e^ h /g 2 is smaller than the critical value r^/v^6, 
there is always a dilatonic black hole solution. The same constraint was derived in 
the case of Einstein-Dilaton-Gauss-Bonnet theory There seems to be, however, 
an additional possibility in the present case. The inequality (47) can be true either 
when Ch > z + > z_ or when Ch < Z- < z + . Since, however, z_ < and Ch > 0, 
the latter case is impossible and only Ch > z + could, perhaps, be realized. Then, 
for x > 1.11949, the inequality < z + < 1 is true. Later we shall justify the 
choice Wh < 1, or equivalently C^ < 1. Thus, solutions could also be present if the 
effective Gauss-Bonnet coupling is in the region " % h > 1.11949 r\. Nevertheless, 
no solutions were found in this new region. 

Before we proceed to describe our numerical procedure and finally plot our solu- 
tions, let us compute analytically the temperature and entropy of the black hole. 
Introducing the Euclidean version of our metric 

ds 2 = e r dr 2 + e A dr 2 + r 2 (d6 2 + sin 2 6d V 2 ) (35) 
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where r is a periodic coordinate ranging over the interval (0,/3), we define the 
temperature as [|12| 



T = — = — - — =^L (36) 

2tt Arc y/g tt g r r \ dr J r=rh An 

where k is the surface gravity of the black hole. On the other hand rearranging the 
equations of motion, we obtain 

d ( e (r-A)/ 2(r _ ^ r 2 _ ^! e (r-A)/ 2[( ^ _ r) (1 _ e - A) + g-A^rV]^ = (37) 
dr V g 2 J 



The equation (|37| ) is an identity. Integrating over the interval (r^,r), we obtain 

2M + D = ^\ 1 (r 2 + —) (38) 

V g 2 

which gives 

T = T SW ^ 2M + V (39) 

\ h g 2 I 

in terms of the temperature of a Schwarzschild black hole T$w — l/87rM = l/Airr^ 
and the asymptotic parameters M and D. Since T is a function of positive constants, 
it never vanishes. As a result, there is no mechanism to prevent the complete evapo- 
ration of the black hole. However, this is indeed the case only when the parameters 
r h and (ph obey the constraint (|34]). If, due to thermal evaporation, the parameters 
of the black hole stop to obey the above constraint, then this gravitational system 
can no longer be described by the concept of a regular black hole. 

The entropy 5* can be derived from the free energy F(f3) of the system as 

S = P&IP-F\ (40) 

The free energy is defined as F = Ie//3, where Ie is the Euclidean version of the 
action 

Ie = J d 4 x^C - J d^x^K (41) 

where K is a suitably subtracted boundary term that comes from R and Rgb- We 
may perform the time and angular integration in ( f4T| ) and by using the equations of 
motion the action takes the form 

I E = 2n(3e^' 2 [ (F - <P') r 2 + ^ (1 - e" A ) (F + 0') - ^ re^r'fi ) |~ (42) 

g g 
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Substituting again the expansions near the horizon as well as near infinity and by 
making use of the definition (HO), we obtain 



4 g 2 r\ 



(43) 



where Ah = ^nr 2 , is the area of the event horizon and S$w = ^-h/^ is the Bekenstein- 



Hawking formula [O for the entropy of the Schwarzschild black hole. 



We note from eq.(|T]) that in the limit <fi — > — oo the effective coupling a'e^/g 2 
vanishes. Then, the contribution of the Gauss-Bonnet term and the SU(2) potential 
to the equations of motion becomes trivial. We can show Q that in this case the only 
acceptable solution is the standard Schwarzschild one with constant dilaton field, 
in agreement with the no-hair theorem. In this limit, we expect the temperature 
T and the entropy S of the coloured black hole to approach the corresponding 
Schwarzschild ones, T S w an d S$w- As we can see from (p9|), since in this limit 
D — > 0, and ([431) , this is exactly what we obtain. 

3. Numerical considerations and conclusions. 



Taking into account the constraints (PB]),(|24]) and ( |3"TD we may conclude that 
the parameters of the problem are r h , <p h , Wh and 71. Note that the equations of 
motion (|6[)-(|lOD do not involve T(r) but only T'(r). The final integration determining 
T(r) involves the integration constant 71 which will fixed by demanding asymptotic 
flatness through fll4]). 



Considering the Yang-Mills equation at the horizon, we obtain 

(e- x )' h w' h + Wh ^ l ~ Wh) = sign w' h = sign w h (w 2 h -l) (44) 
r h 

According to this equation, if we choose the initial value Wh to be greater than one, 
then w' h > 0, which means that w (r) increases with r. If we want to constrain w (r) 
at infinity by w (00) = ±1 then, a local maximum must occur at some point r = r s . 
At this point we would have w (r s ) > 1 and w'(r s ) = 0. Using again the Yang-Mills 
equation, we obtain 

signw"(r s ) = signw (r s ) [w 2 (r s ) — 1] (45) 

If w (r s ) > 1, then w"(r s ) > 0, which means that we can have only a local mini- 
mum at r s . As a result, the initial value of w (r) must fall inside the interval (—1, 1). 
Moreover, we observe that the equations of motion remain invariant under the trans- 
formation w — > — w. Thus, it would be sufficient to choose initial values of Wh in 
the interval (0, 1). 
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The equations of motion (|6"D-(|T0D are invariant under the combined transformation 



0^0 + 0o , r-^r e 0o/2 (46) 

As a result, it is sufficient to vary only one of and <ph. Furthermore, we may 
use the above invariance to set a unique mass scale for all solutions by imposing 
the asymptotic condition cj)^ = 0. This requires a shift <p ^ <p — cf)^ accompanied 
by a rescaling r — > re~^°°^ 2 . Since the radial coordinate has been rescaled, the 
other two asymptotic parameters, M and D, are also rescaled according to the rule 
M — > Me~^°°^ 2 and D — > De^^^l 2 . Similarly, the temperature T is also rescaled as 
T — > Te^°°l 2 . 



In order to perform the numerical integration it is convenient to set % = e = 1. 
We fix the value of the horizon at = 1 and start by giving initial values to 
the remaining parameters <ph and w^. Starting from the expansions (|27| ) - (|30|) , at 
r = r h + e with e ~ O(10 -8 ), we integrate the system (]2Ti|)-(|2~T|) towards r — > oo 
using the fourth order Runge-Kutta method with an automatic step procedure and 
accuracy 10 -8 . The integration stops when the flat space-time asymptotic limit QIBp- 
(|16l) is reached. In the second allowed region of the Gauss-Bonnet effective coupling, 
defined by (47), and for the choices <f)' h = <p' + and (p' h = </>'_ we found solutions with 
regular and non-regular behaviour respectively near the horizon. However, in both 
cases the solutions lacked the right asymptotic behaviour fll3f)-(|i~6l) near infinity. 



Making the choice <p' h = </>'_, which corresponds to the choice (j)' h = <p' + for the 
Einstein-Dilaton-Gauss-Bonnet theory, we were able to find regular asymptotically 
flat black hole solutions in the first allowed region defined by (134]) . For every initial 
value of the shooting parameter 0^ there is a discrete family of initial values of 
Wh which results into a discrete family of asymptotically flat black hole solutions 
characterized by the number n of zeros (nodes) of w(r). We find that as <ph —> — 
oo the variation of the dilaton field <f>(r) with r becomes weaker and the dilaton 
eventually behaves as a constant. The dilaton charge D moves towards zero and the 
horizon takes on its Schwarzschild value = 2M. As we expect, the temperature T 
and the entropy S also take on the corresponding Schwarzschild ones. In the Table 
we display three sets of corresponding values of 4>h, Wh, and (j)^. The displayed values 
of rh, 2M, D, and T are the rescaled ones corresponding to imposing the condition 
0oo = at the end of our computation. Plots involving the dilaton field 0(r), for 
three different coloured black hole solutions, after the imposition of the asymptotic 
condition, are given in Figure 1. Figure 2 depicts the Yang-Mills function w of the 
n = 1,2,3 coloured black hole solutions as well as the metric function gtt{r). The 
dependence of the other metric function g r r{f) is presented in Figure 3. 



S 



TABLE I 



Parameters of solutions for = 1 



4>h 


n 




0oo 


Th 


2M 


D 


T 


-1.0 


1 


0.255948298 


-1.43634 


2.05067 


2.28110 


0.30503 


0.03578 




2 


0.042444284 














3 


0.006916638 












-2.5 


1 


0.262218961 


-2.60493 


3.67836 


3.76249 


0.14371 


0.02123 




2 


0.043590755 














3 


0.007101253 












-5.0 


1 


0.265648181 


-5.00906 


12.2378 


12.2629 


0.04275 


0.00650 




2 


0.044212029 














3 


0.007205569 













In the same region (|34|) and for the choice (p' h = 0', we found solutions which 



had the right asymptotic behaviour (|T^)-(|T^) near infinity and seemed to comprise a 
second group of regular black hole solutions. Examining more closely their behaviour 
near the horizon, they were found to possess a much more complicated structure. 
This structure is characterized mainly by the behaviour of the metric component g rr , 
which is also given in Figure 3, while the other component, g u , exhibits a typical 
black hole behaviour. As we move from infinity towards the origin, the solution 
firstly reaches the value r = r s , around which we may write the expansions 

e- A = X s + X 2 (r-r s ) 2 + ... (47) 

r' = 71 + 72 (r - r s ) + ... (48) 

= <p s + <f)' s {r-r s ) + ^-(r-r s ) 2 + ... (49) 

w" 

w = w s + w' s (r - r s ) + -y-(r - r s ) 2 + ... (50) 

If we substitute these expansions into the equations of motion we obtain a set of 
constraints which determine 72, 4>' s , w' s , 0", and w" as functions of the free parameters 
K, X 2 , 71, 0s, and w s . The curvature invariant R ilvpa R llV(KI near r s was found to be 



= ^ (72 + |) 2 + ^ + / 5)2 4 0( r r, ) ( 5] ) 

z v c r c 



We come to the conclusion that the surface r = r s is a regular one whose existence 
may be interpreted as an unsuccessful attempt of nature to form a regular horizon. 
The significance of this fact will become clear below. Moving beyond the value 
r = r s , the solution eventually stops at the value r = r x . The asymptotic behaviour 
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of the fields near r x is 



e- A = X x + XiVf^+ ... (52) 

T' = 7i + l2V r -r x + ... 

<P = <l) x + (P' x (r-r x ) + ^(r 

w = w x + w' x (r — r x ) + w x (r 

where again 72, <p' x , w' x , (p' x and w x are given through the equations of motion as 
functions of the free parameters \ x , Ai, 71, 4> x and w x . The curvature invariant 
RfivpcrR >lupe7 near r x takes the form 

R^R^ = ( f (72 + ?r^") 2 + T^l — " — <?(^=) - 00 (56) 

V 4 2X X 2r 1 x ) r - r x V r r x 

This means that at the value r = r x the solution ends up in a pure scalar singularity. 
Since no regular horizon exists at some value r = > r x , we may conclude that a 
naked singularity has been formed at r x . 



(53) 

-r s ) 3 / 2 + ... (54) 
~r s f 2 + ... (55) 



For the critical value of (fih, defined by (0), the two solutions and <f>' + coincide, 
which means that the branch of regular black hole solutions meets the branch of the 
solutions that describe a naked singularity. As long as the parameters and (f>h of 
the solution obey the constraint (|34"D, the only acceptable solution of the equations of 
motion is the regular black hole solution (p7|)-(|30D with the flat asymptotic behaviour 
(0)-([16]) near infinity. When and <ph become such - maybe due to thermal 



evaporation - that the constraint ( p4|) is violated , the expansions (27)-(30) break 
down and the solution can no longer be described by the concept of the regular black 
hole. Then, the system shifts to the other branch of solutions, described above, 
corresponding to a naked singularity with exactly the same asymptotic behaviour 
at infinity. 



Alexeyev et al.||, by using a method based on integrating over an additional 
parameter, were able to examine the structure of the black hole solutions found in 
the Einstein-Dilaton-Gauss-Bonnet theory 0] inside the event horizon. According 
to their results, the solution under the regular horizon exist only until the value 
r = r s , where a pure scalar singularity exists. Another solution begins from r s 
which exists only until the "singular horizon" r x . When the Gauss Bonnet effective 
coupling becomes larger, the distance between r s and becomes smaller. Once 
the critical limit (|34j) is reached, = r s = r x and the internal structure vanishes. 
Since the Yang-Mills function w always resembles the behaviour of the dilaton field, 
we expect that the inclusion of the SU(2) Yang-Mills potential respects the above 
internal structure of the regular black hole. Actually, the similarity between the 
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internal structure of the black hole found by Alexeyev et al. and the structure of 
the solution describing a naked singularity found above is obvious. However, the 
interpretation of r x and r s is different: while in the first case they play the role of 
a "singular horizon" and a pure scalar singularity respectively inside the horizon, 
in the second case they stand for a scalar singularity and an unsuccessful event 
horizon respectively. As a result, we can make the further conjecture that as the 
effective coupling becomes larger the internal structure moves towards the horizon, 
for the critical value the points r x , r s and merge, while for values beyond the 
critical point this structure penetrates the horizon and manifests itself as a naked 
singularity. 

Note added in proof. While the present article was being completed, a re- 
lated article titled "Dilatonic Black Holes with Gauss-Bonnet Term" by T. Torii, 



H. Yajima and K. Maeda, |gr-qc/9606034| came into our attention. 



Acknowledgements. We wish to thank J. Rizos and N. Mavromatos for many 
helpful discussions. 
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Figure 1: Dilaton field for coloured black hole solutions. Each curve corresponds to 
a different solution characterized by a different initial value of 4>h- 
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Figure 3: Dependence of the metric component g rr for r^ — l and 4>h — ~ 2- While 
for <j)' h = g rr moves towards an infinite value leading to the formation of an event 
horizon, for <p' h = <p' + it reaches a maximum value at r s and then declines leaving 
unshielded a scalar singularity. 
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